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DRAG OF A SUPERCRITICAL BODY OF REVOLUTION IN FREE FLIGHT AT 
TRANSONIC SPEEDS AND COMPARISON WITH WIND-TUNNEL DATA 


By J. W. Usry and John W. Wallace 
Langley Research Center 

SUMMARY 

The forebody drag of a supercritical body of revolution was measured in free flight 
over a Mach number range of 0.85 to 1.05 and a Reynolds number range of 11.5 x 10° to 
19.4 x 10 6 and was compared with wind-tunnel data. The forebody drag coefficient for a 
Mach number less than 0.96 was 0.111 compared with the wind-tunnel value of 0.103. A 
gradual increase in the drag occurred in the Langley 8 -foot transonic pressure tunnel at a 
lower Mach number than in the Langley 16 -foot transonic tunnel or in the free-flight test. 
The sharp drag rise occurred near Mach 0.98 in free flight whereas the rise occurred near 
Mach 0.99 in the Langley 16 -foot transonic tunnel. The sharp rise was not as pronounced 
in the Langley 8-foot transonic pressure tunnel and was probably affected by tunnel -wall- 
interference effects. The increase occurred more slowly and at a higher Mach number. 
These results indicate that the drag measurements made in the wind tunnels near Mach 1 
were significantly affected by the relative size of the model and wind tunnel. 

INTRODUCTION 

A free-flight test was made to provide data at high Reynolds numbers on an aero- 
dynamic configuration designed to fly in the sonic or near -sonic speed regime. The flight 
test provided sting -free, wall -interference -free data on the transonic drag characteristics 
of a low -drag body of revolution. These data along with transonic wind-tunnel data mea- 
sured at the same Mach numbers and Reynolds numbers are being used as part of the data 
base in applications such as advanced technology aircraft. One approach to the design of 
these types of aircraft is to define the aerodynamic configuration that allows the aircraft 
to cruise near sonic speeds without large penalties in performance due to the drag-rise 
characteristics. This type of design would require a low subsonic drag configuration with 
a drag divergence Mach number near 1. 

Wind-tunnel studies in the Langley 8-foot transonic pressure tunnel (TPT) have pro- 
duced several bodies of revolution with low drag which appear to have the desired drag 
divergence characteristics. Acquisition of accurate reliable data near Mach 1 is difficult, 
however, because of effects due to the tunnel boundaries. For this reason, a free-flight 


test of one configuration tested in the Langley 8 -foot transonic pressure tunnel and in the 
Langley 16-foot transonic tunnel (TT) was made for data correlation purposes. 

The purpose of this paper is to present the results of the flight test and to compare 
the data with wind-tunnel results. Specifically, forebody drag and base pressure coeffi- 
cients are presented and compared with wind-tunnel data over a Mach number range from 
0.85 to 1.05 and a Reynolds number range from 11.5 x 10® to 19.4 x 10®. A comparison 
of the drag-rise trends near Mach 1 determined from free-flight and wind-tunnel data is 
presented and the effects due to test-environment differences are considered. Also his- 
tories of the flight-test environment parameters and onboard measurements are pre- 
sented. The test was conducted at the NASA Wallops Station, Wallops Island, Virginia. 

SYMBOLS 


a x 

c A,b 

c A,fb 

A CA,fb 

C A)t 

C D,f 

c p,b 

c 

d b 

d max 

g 

h 


accelerometer reading along body X-axis, g units 

Su 

base axial drag coefficient, -C« b tt 

forebody axial -force coefficient, t - Ca b 

incremental forebody axial -force coefficient, CA ? fb - (CA,fb)jy[_o 9 

Wa 

total axial -force coefficient of model, £ 

’ q S 

^OO 

skin-friction drag coefficient, based on S 

p. - p 

base pressure coefficient, - 2 - — — 

” OO 

chord, centimeters 

diameter at base of model, centimeters 
maximum diameter of model, centimeters 
acceleration due to gravity, m/sec^ 
altitude, meters 
airfoil thickness, centimeters 
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h r relative humidity, percent 

l length of model, centimeters 

M Mach number 

o V l 

Np_ Reynolds number, °° °° 

M- oo 

base pressure, N/m2 
p c saturation pressure, N/m^ 

P t total pressure, N/m^ 

P v vapor pressure, N/m2 

p^ free-stream static pressure, N/m^ 

Ap differential pressure, N/m^ 

Ap^ differential pressure in pitch plane, N/m^ 

Ap^ differential pressure in yaw plane, N/m2 

free-stream dynamic pressure, N/m^ 

R gas constant for air 

r body radius, centimeters 

S reference area based on maximum diameter of model, meters^ 

reference area based on diameter at base of model, tfdb/4, meters^ 

T temperature, K 

t time, seconds 

V earth relative velocity, m/sec 
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local speed of sound, m/sec 


V s 

V w wind velocity, m/sec 

Voo free-stream velocity, m/sec 

W weight of model, 36.35 kgf or 356.47 newtons 

X,Y,Z body-axis system 

x,y 

x i 
a 

fi 

V 

e 

Moo 

Poo 
<7 

i \f/' = xfs - ^ w , degrees 

\p w azimuth from true north of wind, degrees 

RESEARCH MODEL DESCRIPTION 

A sketch of the model with dimensions and other details is presented as figure 1 and 
photographs of the model on the aircraft are shown as figure 2. The aircraft support rig 
was designed to carry the model so that no protrusions or indentations appeared on the 
surface of the model. The model was a smooth aerodynamic body of revolution. The 


model coordinates, centimeters 

any variable listed in table III 

angle of attack, degrees 

angle of sideslip, degrees 

, 9 oN 1 / 2 

resultant angle, [a* + , degrees 

elevation angle, degrees 
free-stream viscosity, N-sec/m 2 
free-stream density, kg/m 3 

standard deviation of variable denoted as subscript 
azimuth from true north of payload, degrees 
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exterior shape was defined to provide low drag characteristics with drag divergence near 
M = 1. The body radius as a function of longitudinal body station is presented as table I. 

The length of the model was 114.30 centimeters and the maximum diameter was 
12.70 centimeters. The diameter of the model at the base was 3.28 centimeters. The 
model was aerodynamically stabilized during the flight by using fins which had a biconvex 
airfoil section with a thickness -chord ratio of 0.03. The fins were swept 45° and the root 
chord projected to the body center line was 17.20 centimeters in length. The area rule 
was used to provide a smooth distribution of the longitudinal cross-sectional area in the 
vicinity of the fins. Fin dimensions and a set of typical airfoil coordinates are listed in 
table n. The center of gravity was located 53.46 centimeters aft of the nose tip (46.7 per- 
cent of the model length). Differential-pressure ports in the pitch and yaw planes were 
located on the model as shown in figure 1. A total-pressure port was located at the 
nose tip and a base -pressure port was located midway between two fins at a distance of 
1.32 centimeters from the model center line. A transition strip was located 2.54 centi- 
meters aft of the nose. The strip was 0. 127 centimeter wide in the stream direction and 
consisted of No. 120 carborundum grit particles in a concentration of approximately 
20 per centimeter. Transition strip location and grit characteristics were the same 
as those used on the wind-tunnel model and were determined by using the method of 
reference 1. 


RESEARCH MODEL INSTRUMENTATION 

An FM/FM telemeter transmitted six channels of information to ground receiving 
stations. These transmissions included total pressure, two differential pressure mea- 
surements on the nose, two longitudinal accelerations, and base pressure. The following 
table lists the parameters measured and the range of each instrument used. The expected 
accuracy (la deviation) of each measurement after data reduction to engineering units was 
assumed to be ±1 percent of the range of the instrument. An error analysis is presented 
in the appendix. The orientation of the body-axis system for the results presented in this 
paper is shown in figure 3. 


Measurement 

Range of 
instrument 

Total pressure, kN/m^ 

Differential pressure (two), kN/m2 

Longitudinal acceleration, g units 

Base pressure, kN/m2 

0 to 172.4 
0 to 68.9 
f0.25 to -1.00 
\0.00 to -0.50 
0 to 103.4 


5 



TEST ENVIRONMENT 


Trajectory Parameters 

The test environment was achieved by dropping the research model from an aircraft 
at an altitude of 9411 meters and a velocity of 182.9 m/sec. The mission profile is shown 
in figure 4 as a plot of altitude against horizontal range. The time at release from the 
aircraft was taken as t = 0. Time after release and the Mach number range of interest 
are indicated on the plot. These data were calculated from the best available radar track 
of the model. 

Histories of the altitude, velocity, flight -path angle, Mach number, dynamic pres- 
sure, and Reynolds number are shown in figure 5 for the prime data period. Altitude and 
flight-path angle were calculated from FPQ-6 radar positional data. Free-stream veloc- 
ity was obtained from FPQ-6 radar positional data, FPQ-6 radar Doppler data, and 
Jimsphere wind data. Mach number, dynamic pressure, and Reynolds number were cal- 
culated from free-stream velocity, atmospheric speed of sound, temperature, and density. 
Mach number was also calculated by using the tables of reference 2 and the ratio of the 
total and atmospheric pressure. Accuracies of these parameters and the method used to 
estimate the accuracy are presented in the appendix. 

Atmospheric Parameters 

Variations of the atmospheric temperature, pressure, density, speed of sound, and 
relative humidity with altitude are presented in figure 6. The temperature, pressure, and 
humidity were measured by use of a standard radiosonde balloon launched at the same 
time that the model was released from the aircraft. The measurements, therefore, were 
made within 25 minutes after release. The other parameters were calculated from these 
measurements. In addition, the parameters from the 1966 U.S. Standard Atmosphere 
Supplements (ref. 3) are shown for comparison. Accuracies of the measured quantities 
are based on those in reference 4. 

Variations of the atmospheric wind velocity and direction with altitude are shown in 
figure 7. These data were measured by using the FPS-16 radar -Jimsphere balloon sys- 
tem and have a root-mean-square accuracy of 0.5 m/sec (see ref. 5). The Jimsphere 
balloon was launched about 25 minutes after the model was released from the aircraft. 
Therefore, these data were obtained within 50 minutes after model release. These data 
were used in the calculation of free-stream velocity, Mach number, and dynamic pressure. 

RESULTS AND DISCUSSION 

Histories of the longitudinal acceleration, total and base pressures, and the differ- 
ential pressures measured onboard the model are shown in figure 8. The measurements 
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made by the two longitudinal accelerometers were essentially the same and the curve 
in figure 8(a) represents both measurements. The primary data period was from 
t = 22 seconds to t = 36 seconds corresponding to M = 0.85 to 1.05 and Reynolds 
numbers from 11.5 x 10® to 19.4 x 10®. 

Angle -of -Attack Effects 

A model geometrically similar to the flight model was tested in the Langley 8 -foot 
transonic pressure tunnel and the Langley 16 -foot transonic tunnel to determine the vari- 
ation of the differential pressures with angle of attack and Mach number. These varia- 
tions were linear with angle of attack for a given Mach number for the range of angles 
and Mach numbers of interest. The slopes of these curves (differential pressure coeffi- 
cient per degree) plotted against Mach number are shown in figure 9. The difference in 
3^Ap/q oo)/ 9 ** between the two sets of wind-tunnel data at M = 1 is equivalent to a dif- 
ference in angle of attack of 0.05° for the maximum angles obtained in flight; thus, the 
faired curve of figure 9 was used with the differential-pressure data of figure 8 measured 
in flight to estimate the angles of attack and sideslip. 

Time histories of the angles of attack, sideslip, and resultant angle are shown in 
figure 10. These data show that the model flew at a trim angle of attack which decreased 
with time or increasing Mach number. As the model accelerated through M = 1, the 
trim angle increased to 2° but continued the same decreasing trend after passing through 
M = 1.0. (The estimated standard deviation of the resultant angle at M = 1 was 0.4°.) 
The reason for the model flying at these trim angles is not apparent and cannot be deter- 
mined from the limited measurements made. 

The variation of axial-force coefficient with angle of attack was measured in the 
wind tunnels at an angle of sideslip of 0° and is shown in figure 11. These data were used 
to obtain the variation of axial-force coefficient with Mach number at an angle of attack 
of 0° for comparison with the free -flight test results. In addition, these data are pre- 
sented to show trends of with angle of attack and to illustrate that 

decreases with increasing angle of attack for small angular deviations. 

Comparison of Flight Data With Wind-Tunnel Data 

The forebody axial -force and base pressure coefficients at the flight trim angles 
are compared with wind-tunnel data at an angle of attack of 0° in figures 12 and 13. The 
flight data were calculated as follows: 

C A,fb = C A,t - (Base axial drag coefficient) 
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where 


C P,b 


- P 


00 


q 


oo 


and 


C A,b “ _c p,b y 


The 16 -foot -tunnel data in figures 12 and 13 were measured at Reynolds numbers 
from 15.6 x 10® at M = 0.85 to 17.04 x 10® at M = 1.01 and the 8-foot-tunnel data were 
measured at a Reynolds number of 16 x 10® for all Mach numbers. The flight Reynolds 
numbers varied from 11.5 x 10® at M = 0.85 to 19.4 x 10® at M = 1.05. (See fig. 5.) 

The from free -flight data (see fig. 12) for M < 0.96 was 0.111 whereas 

the value measured in both wind tunnels was 0.103. The error analysis presented in the 
appendix shows that one-half of this difference could be attributed to the estimated error 
in C A ffo. The Reynolds number varied during the flight but was essentially constant for 
both wind-tunnel tests. Also, sting -interference effects could have caused the pressure 
distribution on the wind-tunnel model afterbody to be considerably different from that on 
the flight model. Both of these effects could have contributed to the drag level differences 
for M < 0.96. 

Figure 13 shows the base pressure coefficient as a function of Mach number with 
the wind-tunnel data shown for comparison. The base pressures were higher in both tun- 
nels for M < 0.96 and indicated that the pressures on the afterbody of the wind-tunnel 
model could also have been higher. This condition would make the wind-tunnel 
lower than the free -flight values since these pressures were positive. These differences 
in base pressure alone, however, represent 2.4 percent of the total drag at M = 0.85 and 
become practically insignificant above M = 0.96. The sting configurations used in the 
wind-tunnel tests, however, were designed to minimize sting interference effects by use 
of the method of reference 6. 

Differences which can be attributed to Reynolds number effects are indicated by fig- 
ure 14 where a theoretical estimate of the skin-friction drag coefficient is shown as a 
function of Mach number. These coefficients were calculated by using the Sommer -Short 
T' method and the Karman-Schoenherr flow equation for a turbulent boundary layer and 
the Blasius- Sutherland equation for a laminar boundary layer. (See refs. 7 and 8.) For 
this analysis the flow was assumed to be laminar forward of the transition strip and turbu- 
lent aft of the strip. Flight conditions were used to obtain the solid curve in figure 14 
whereas the broken curve was obtained by using a constant Reynolds number and stagnation 
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temperature for all Mach numbers and thus simulates the wind-tunnel conditions. Both 
curves at the top of the figure were calculated by assuming turbulent flow on the body 
(aft of transition strip) including the fins. The broken curve in the bottom of the figure 
was obtained by using flight conditions and laminar flow over the fins. These curves 
indicate that the drag coefficient obtained in flight should be slightly higher than the wind- 
tunnel values for M < 0.96 and slightly less for M > 0.97, provided the flow was turbu- 
lent over the fins. The actual flow probably was mixed since local Reynolds numbers 
based on fin chord varied from about 1 x 10® at the fin-body juncture to 0 at the fin tip at 
the lower Mach numbers and from about 2 x 10® to 0 for the larger Mach numbers. Also, 
parts of the fins were immersed in the turbulent boundary layer over the body. The 
model tested in the wind tunnels was geometrically similar and should have had the same 
flow characteristics as the flight model had, provided the boundary and sting interference 
effects did not affect the nature of the flow. These data do indicate, however, that con- 
siderable differences between the flight and wind-tunnel data could be attributed to 
Reynolds number and flow character effects. 

A comparison of the incremental axial-force coefficients is shown in figure 15 as 
ACA,fb plotted against Mach number where AC^f^ is the difference in C^ft, from 
the value at M = 0.9. The gradual increase in drag obtained in each of the three tests 
is probably due to the initial formation of weak shocks on the aft end of the body and the 
resulting upstream pressure adjustment through the subsonic part of the boundary layer. 
This gradual increase began at a lower Mach number in the 8 -foot transonic pressure tun- 
nel than in either the free -flight test or the 16 -foot transonic tunnel and was probably due 
to boundary interference effects. The AC^fb data from the free flight and 16 -foot 
transonic tunnel agree very well up to M = 0.98. At this point the free -flight drag data 
increased sharply whereas the sharp rise in the 16 -foot -tunnel data occurred near 
M = 0.99. The sharp rise in the data from the 8-foot transonic pressure tunnel was not 
as pronounced. The increase in drag occurred more slowly and at a higher Mach num- 
ber. The generally good agreement between the flight data and the 16 -foot -tunnel data 
and the fact that the 16-foot- and 8-foot -tunnel data were obtained under very nearly iden- 
tical conditions (except for tunnel size) indicate that the drag data obtained near Mach 1 
on this model may be significantly affected by factors related to the relative size of the 
model and the tunnel. 


CONCLUDING REMARKS 

Drag characteristics of a body of revolution with cruciform fins were determined 
in free flight over a Mach number range of 0.85 to 1.05 and a Reynolds number range of 
11.5 x 10® to 19.4 x 10® and were compared with data measured in the Langley 8-foot 
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transonic pressure tunnel and the Langley 16 -foot transonic tunnel. The free -flight 
model was dropped from an altitude of 9411 meters at an initial velocity of 182.9 m/sec. 

Differences in the flight and wind-tunnel data could be attributed to the combined or 
individual effects of Reynolds number differences between the flight test and the wind- 
tunnel tests, flow quality conditions, boundary effects, and instrumentation differences. 
Further wind-tunnel studies and experimental tests are needed to determine the effects 
of each of these parameters on the drag. 

The axial -force coefficient for a Mach number less than 0.96 was 0.111 in free 
flight compared with 0. 103 for the wind tunnel. A gradual increasing trend in the drag 
occurred in the Langley 8-foot transonic pressure tunnel at a lower Mach number than 
in the Langley 16 -foot transonic tunnel or in the free -flight test. The sharp drag rise 
occurred near Mach 0.98 in free-flight data and near Mach 0.99 in the Langley 16-foot 
transonic tunnel data. The sharp rise was not as pronounced in the data from the Langley 
8-foot transonic pressure tunnel. The increase in drag occurred more slowly and at a 
higher Mach number. These results indicate that the drag measurements made in the 
wind tunnels near Mach 1 were significantly affected by the relative size of the model 
and wind tunnel. 

Langley Research Center, 

National Aeronautics and Space Administration 
Hampton, Va., October 27, 1971. 
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APPENDIX 


ERROR ANALYSIS OF FREE -FLIGHT DATA 

The root -sum -square (RSS) method for combining errors from independent sources 
was used to estimate accuracies of the free -flight data. The major parameters of inter- 
est are the forebody axial-force coefficient C^fb* Mach number M, and Reynolds num- 
ber Npjg. The assumption was made that the individual error source contributions to 
the parameters are statistically independent and normally distributed with mean zero. 
The RSS method may be expressed as follows: 



where 

CTy one standard deviation in y = f^XjJ 

CT Xi one standard deviation in independent variables Xj (i = 1, 2, 3, . . ., n) 

partial derivative of y = ffa) with respect to x< 

oX| ' ' 

The study performed here is concerned with the accuracy of axial -force coefficient, Mach 
number, and Reynolds number. 

The forebody axial -force coefficient is a function of several independent variables 
and may be expressed as follows: 


c A,fb = -■ 


2RgT[wa x -(p b -p oo )S b ] 


s|p oo -18.147h r |2.685 + 3.537xl0" 3 ||(T-459.67)J 2 ’ 245 ^[v+V w (cos i//-^ w )cos 0 ] 


(A2) 

By applying the RSS method (eq. (Al)) to equation (A2), the standard deviation of C ^ ^ 
can be estimated as follows: 


CT CA,fb ” 




(A3) 
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APPENDIX - Concluded 


Histories of the standard deviations were calculated by using the chain-rule method 
illustrated in figure 16. Values of the standard deviations used in these calculations are 
listed in table HI. The expected accuracy (lcr deviation) of each onboard measurement 
after data reduction to engineering units was assumed to be ±1 percent of the range of the 
instrument. 

Standard deviations of the wind data ^cry and cr^ j and atmospheric data (oj,, 

Ok . and o n ^ were obtained from references 4 and 5. Standard deviations for radar 

n r ^coy 

data (cry, cr^, and o^j were obtained from a radar study conducted at Wallops Island, 
Virginia. 

Table IV presents the values and the estimated standard deviations of M, 

and Nj^g at M = 1. The values and estimated standard deviations for the intermediate 
variables in these estimations are also listed. The standard deviation of C^fb is 
about 2 percent of the calculated value at M = 1.00. The standard deviation is less than 
0.2 percent for M and less than 0.4 percent for N^g. 

Figure 17 shows the variation of the error contribution of four of the basic vari- 
ables and On with M. The four variables make the largest error contributions 
~A,fb 

to On and the error contribution of the a x measurement is an order of magnitude 

^A,fb 

larger than the other three error sources and clearly illustrates the need to improve the 
accuracy of the complete data acquisition system. 

Figure 18 shows the variation of ^ and ± with M. Also 

indicated in the figure is M ± at M = 0.97. The error in at M * 0.80 

is 8 percent of and decreases with increasing M. At M = 1, the error is 

2 percent. 
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TABLE II.- FIN DIMENSIONS 


The coordinates of any airfoil section in the chordwise direction are given by 




where -g= = 0.03. 


c = 12.70- 


/ V.50 

L 



c = 17.20 


Typical coordinates of airfoil section, c = 12.70 cm 


x, cm 

y, cm 

x, cm 

y, cm 

x, cm 

y, cm 

x, cm 

y, cm 

0.0000 

0.0000 

3.3020 

0.1466 

6.6040 

0.1902 

9.9060 

0.1308 

.2540 

.0149 

3.5560 

.1537 

6.8580 

.1892 

10.1600 

.1219 

.5080 

.0292 

3.8100 

.1600 

7.1120 

.1877 

10.4140 

.1125 

.7620 

.0429 

4.0640 

.1659 

7.3660 

.1857 

10.6680 

.1024 

1.0160 

.0561 

4.3180 

.1709 

7.6200 

.1829 

10.9220 

.0917 

1.2700 

.0686 

4.5720 

.1755 

7.8740 

.1796 

11.1760 

.0805 

1.5240 

.0805 

4.8260 

.1796 

8.1280 

.1755 

11.4300 

.0686 

1.7780 

.0917 

5.0800 

.1829 

8.3820 

.1709 

11.6840 

.0561 

2.0320 

.1024 

5.3340 

.1857 

8.6360 

.1659 

11.9380 

.0429 

2.2860 

.1125 

5.5880 

.1877 

8.8900 

.1600 

12.1920 

.0292 

2.5400 

.1219 

5.8420 

.1892 

9.1440 

.1537 

12.4460 

.0149 

2.7940 

.1308 

6.0960 

.1902 

9.3980 

.1466 

12.7000 

.0000 

3.0480 

.1389 

6.3500 

.1905 

9.6520 

.1389 






TABLE m.- ASSUMED VALUES OF STANDARD DEVIATIONS (la) FOR CALCULATIONS 


<7ax> S °- 005 

ap^, kN/m^ 1.034 

a n , kN/m2 1.724 

p t 

a v > m/sec 0.5000 

v w 

CT »/V rad 0.01745 

<jy, m/sec 0.3048 

rad 0.00015 

c tq , rad 0.00015 

a T , K 0.70 

% for ~ 

273.15 K > T = 233.15 K 0.033 

T 5 273.15 K 0.017 

a p for - 

h < 3048 m 0.023 

3048 m < h < 6096 m 0.033 

6096 m < h < 9144 m 0.040 
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Photographs of model on aircraft. 




(b) Rear quarter view. 
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Horizontal range, k* 

Figure 4.- Variation of altitude with horizontal range based on radar data. 
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Time, sec 


(a) Time histories of altitude, velocity, and flight-path angle. 
Figure 5.- Time histories of test conditions. 
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Temperature, °K 


(a) Temperature. 


Pressure, kfi/m 


(b) Pressure. 
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(c) Density. (d) Speed of sound. 
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Relative humidity, percent 


(e) Relative humidity. 

Figure 6.- Atmospheric properties as a function of altitude. 








(a) Longitudinal acceleration. 
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(b) Total and base pressures. 
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(c) Differential pressures in pitch and yaw planes. 
Figure 8. - Time histories of data from onboard instruments. 
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Figure 11.- Variation of forebody axial -force coefficient with angle of attack from the 
Langley 8 -foot transonic pressure tunnel and the Langley 16 -foot transonic tunnel. 
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Figure 12.- Variation of forebody axial -force and base drag coefficients 

with Mach number. 
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Figure 15.- Comparison of incremental axial -force trends for data 
from tests in the Langley 8 -foot transonic pressure tunnel, the 
Langley 16 -foot transonic tunnel, and in free flight. 
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m of drag coefficient and standard deviation with Mach number. 
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